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Abstract: Modified nucleic acid bases are most commonly found in tRNA. These may contain modifications from simple 
methylation to addition of bulky groups* Methylation of the four canonical nucleotide bases at a wide variety of positions is 
particularly prominent among the known modification* Methylation of N2 group of guanine is a relatively common modification 
in tRNA and rRNA. N^-methylguanosine (m^G) is the second most often encountered nucleoside in coli tRNAs* N^, N^- 
dimethylguanosine (m^zG) is found in the majority of eukaryotic tRNAs and involved in forming base pair interactions with 
adjacent bases. Hence, in order to understand the structural significance of these methylated nucleic acid bases we have carried out 
molecular dynamics simulation to see the salvation effect. The results obtained shows iso-energetic conformational behaviors for 
m^G and m^zG, The simulation trajectory of m^G shows regular periodical fluctuations suggesting that m^G is equally stable as 
either s-cis or s-trans rotamers. The two rotamers of m^G may interact canonically or non-canonically with opposite base as s-trans 
m^G26:C/A/U44 and s-cis m^G26:A/U44, The free rotations around the C-N bond could be the possible reason for these iso- 
energetic conformations, Dimethylation of G has almost no influence on base pairing with either A or U, Thus, these results reveal 
that modified nucleosides m^G and m^zG may play an important role to prevent tRNA from adopting the unusual mitochondrial 
like conformation. 
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Introduction 

RNA molecules undergo extensive post-transcriptional 
modifications that are important for their biological activities. Post- 
transcriptional modifications have been known as a natural 
mechanism to provide structural stability across the wide range of 
temperature in archaea as well as bacteria [I]. Transfer RNAs have 
the largest number and the greatest diversity of modifications: base or 
ribose methylation, base isomerization, base reduction, base thiolation 
and more complex hypermodifications [2, 3]. An important 
characteristic of tRNA is the presence of high content of modified 
nucleosides of which methylation represents the principle post- 
transcriptional modification during its maturation. In the maturation 
process of tRNA, transfer of methyl group occurs at polynucleotide 
level through an S-adenosyl-L-methionine donor, resulting in 
modification of heterocyclic base, the ribose moiety, or both [4]. The 
family of structurally related nucleosides m^G, m^zG, m^Gm and 
m^zGm, are from known archaeal tRNA sequences. These modified 
nucleosides are conserved at only two locations, position 10 first base 
in the proximal position of the dihydrouridine (DHU) arm and at 
position 26, junction between the D-stem and the anticodon stem, 
where they play crucial roles in the control and stabilization of the 
tertiary L fold structure of the tRNA [5,6]. The m^G and mhG 
modifications in tRNA are found not only at position 26 but also at 
positions 6, 7,9, 10, 18 and 27 in various organisms [7]. 
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Experimentally, it has been found that level of certain modified 
nucleosides in archaeal thermophiles play major stabilizing role 
beyond the effects of magnesium ion binding and G-C content of 
tRNA [8]. Earlier study involving three-dimensional models of yeast 
tRNA^^' derived from X-ray crystallographic data implies that 
m^2G26 functions as a molecular hinge. 

This hinge adjusts the angular position of the D-stem and the 
anticodon stem during protein synthesis, thus maintaining a certain 
rigidity/ flexibility in this part of tRNA [9]. Nuclear magnetic 
resonance studies on the resonance of the methyl proton in yeast 
tRNA^^' also provide evidence to support the notation that m^2G26 
has a significant role in regulating the stacking and conformational 
dynamics of this region of tRNA molecule [10]. The yeast tRNA 
(m^2G26) methyl transferase is dependent on the D-stem sequence and 
size of variable loop for the synthesis of N^-N^ dimethyl guanosine at 
26^^ position [II]. Mutations were introduced in both the D-stem 
and the variable loop of tRNA^'P to obtain dimethylation of the 
normally unmodified G26 by the yeast N^, N^-dimethyl G26- 
methyltransferase [12]. The presence of m^2G26 in cytosolic tRNA 
may avert the molecule from adopting an unusual mitochondrial 
tRNA pattern folding and instead, allow it to fold into the canonical 
cloverleaf model. Through screening of the tRNA sequence and gene 
database it was revealed that some cytosolic tRNAs have the potential 
to fold into alternate structures. It was further noted that when a 
tRNA had the potential for this alternate folding, m^2G was found at 
position 10 and 26 presumably to block the formation of this non- 
standard folding pattern [13]. The methylated guanosine from 26 
position of tRNA may have role in regulating the stacking 
interactions and the conformational dynamics [14]. 

N^-methylguanosine is found in both helical and loop regions of 
RNA secondary structure [15, 16] and it can exist in either s-cis or s- 
tmns rotamers [18]. Incorporation of m^G was found to be iso- 
energetic with G in the duplex context as well as in GNRA (N = any 
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nucleotide and R = purines) tetra loops [17]. The two ro tamers of 
m^G, s-cis and s-trans have been found equally stable in RNA duplex 
[17] and in tRNA [18]. Free rotations around C-N bond aids m^G to 
get energetically two stable conformations and to form base pair 
interactions such as s-cis m^G:A/U or s-trans m^G:C/A/U [18]. 
This phenomenon has not been studied in detail at atomic level. 
Hence, present attempt has been made to investigate the dynamic 
behavior of iso-energetic conformations of m^G/m^zG using MD 
simulation technique. The results clearly show that the iso-energetic 
nature of s-trans and s-cis conformations could be because of free 
rotations around the C-N bond. Before this various computational 
techniques have been used to understand the conformational behavior 
and dynamics of many complex modified nucleosides [18-23]. 
However, all-atom molecular dynamics simulations of entire solvated 
ribosome, mRNA and tRNA complex have been studied to find out 
motion of tRNA from the A/T state into the A site [24]. 
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Figure 1. The cloverleaf structure showing m^G and m^2G at 26^^ position. 



Computational details 

Molecular dynamics (MD) simulations were performed using 
Amber 10 simulation suite on (HP ProLiant-MLI50G6) server in 
order to highlight the influence of explicit solvation on the 
conformation of modified nucleosides N^-methylguanosine (m^G) 
and N^-N^ dimethyl guanosine (m^zG). PCILO predicted preferred 
conformations of m^G and m^zG [18] were used as starting 
geometries for MD simulation studies. Antechamber software was 
used to calculate charges. Each nucleoside was solvated by 729 
SPC/E water molecules fillmg a 34.20 X 31.05 X 31.85 
rectilinear box with I.O water density [25]. Simulations were 
performed under periodic boundary conditions by employing the 
Particle Mesh Ewald [26] method to calculate long-range 
electrostatic interactions. MD trajectories were propagated at 2.0 fs 
time step using the shake algorithm [27] to all hydrogen atoms with 



non-bonded cutoff of 10 A. The non-bonded pair list was updated 
by every 10 steps. The trajectories were calculated by maintaining 
constant temperature (300 K) and constant pressure (latm) at 2 fs 
time step according to Berendsen coupling algorithm [28]. 

An equilibration protocol similar to the earlier molecular 
dynamics simulation study of nucleic acids was applied [29, 30]. The 
equilibration protocol consisted of 5000 steps of steepest descent 
minimization followed by 5 ps of MD at 300 K applied to relaxation 
of initial strain present between water molecules and N^-methyl 
derivatives of guanosine. In the next step N^-methyl derivatives were 
fixed while water molecules were allowed to relax at 100 K (I ps), 
200 K (I ps), and finally at 300 K for 198 ps, thus equilibration 
protocol was completed at 200 ps. 

Equilibrated system was further subjected to 5000 steps of 
steepest descent minimization to remove bad contacts between water 
molecules and nucleic acid bases. In further steps of MD simulation, 
no positional constrains were applied to the system and the 
temperature was progressively increased to 300 K in steps of 50 K 
with I ps at each step. Finally system was subjected to production 
MD of 20 ns at 300 K temperature and constant pressure (latm) 
with fully solvated and neutralized system. PTRAJ module of Amber 
Tool 10 was used for analysis of average structures [31]. 

Results and Discussion 

Dynamic behavior of JSF -methylguanosine ( m^G) 

The preferred conformation of N^-methylguanosine [18] (Fig. 
2A) has been used as a starting geometry for 20 ns molecular 
dynamics simulation study. In order to confirm the iso-energetic 
conformational behavior of m^G we have analyzed four different 
average structures particularly at 0 to I ns (Fig. 2B), 3 to 4.5 ns (Fig. 
2C), 5 to II ns (Fig. 3B) and 13 toI9 ns (Fig. 3C) and three 
snapshot structures particularly at 2ns (Fig. 3 A), 1 2ns (Fig. 2D) and 
20ns (Fig. 2E) of 20ns total simulation period. The geometrical 
parameters are mentioned in table I. The selection of average and 
snapshot structures have been made based on the conformational 
flexibility observed during the MD simulation trajectory (Fig 4, 6) 
similarly as per our earlier conformational studies of yW [20], 
OHyW [21] and ac^C [32]. 

Stabilization of s-trans m^G26 conformation 

The MD simulation average structures for m^G taken at 0 to I ns 
(Fig. 2B), 3 to 4.5 ns (Fig. 2C), and snapshot structures at 12 ns (Fig. 
2D) and 20 ns (Fig. 2E) shows the "proximal'' or s-trans orientation 
with imidazole ring of guanosine. 

The Methyl group of m^G point towards the N(3) atom of 
guanosine as observed in earlier study [18]. This s-trans or proximal 
orientation would allow Watson-Crick base pairing of m^G26 with 
C44 and non Watson-Crick base pairing with A/U44 at the hinge 
region of tRNA. Similar kind of s-trans orientation for m^G has been 
observed in our earlier conformational study [18] along with crystal 
conformer of m^GIO [33], where it forms Watson-Crick base pairing 
with C25. 

The average structure obtained at 0 to I ns (Fig. 2B) maintains 
the initial geometry (Fig. 2A) [18] by preserving s-trans or 
''proximaF' conformation for N^-methyl substituent of guanosine 
(26^^), which is stabilized by hydrogen bonding interaction between 
N(3)...HC(I0) (Fig. 2B and Table I). This average structure (Fig. 
2B Table I) shows deviations for torsion angle p by 87° and % by 
155° whereas a retains its initial geometry [18] as found in crystal 
conformer lEHZ.pdb [33] and 6TNA.pdb [9]. A large deviation 
around the torsion angle p is due to rotations around C-N bond. 
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Figure 2. (A) PCILO predicted most stable structure of N -methylguanosine (m G) [18]. The methyl group has been given 'Yellow' colour for clear identification 
(B) Average structure of N^-methylguanosine for 0-1 ns. (C) Average structure of N^-methylguanosine for 3-4.5 ns. (D) Snapshot structure of N^- 
methylguanosine at 12 ns. (E) Snapshot structure of N^-methylguanosine at 20 ns. 



Next average structure taken at 3 to 4.5 ns (Fig. 2C) is also stabilized 
by hydrogen bonding between N(3). . .HC(IO) along with this, 
interaction between N(3)...HC2' (Fig. 2C and Table I) provides an 
additional structural stability to this average structure (Fig. 2C), as 
observed in our earlier conformational study of m^G [18]. The 
torsion angle p shows large deviation (120°) whereas a maintains 
starting value as compared with initial structure of m^G [18]. 

Snapshot structure selected at 12 ns (Fig. 2D) prefers s-trans or 
^proximaF' conformation for m^G and stabilized by N(3). . .HC(IO), 
N(3)...HC2' and N(3)...H02' mteractions (Fig. 2D and Table I) 
similar to earlier results of m^G [18]. This snapshot structure shows 
similar conformation for torsion angle a while torsion angles p and % 
deviates to large extent from initial structure (Fig.2A) as observed in 
crystal structure lOBS.pdb [34]. 

Second snapshot structure (Fig. 2E and Table I) selected at final 
trajectory (20 ns) of simulation study also preserves s-tmns 
conformation for m^G and stabilized by intramolecular interactions 
between N(3). . .HC(IO), N(3)...HC2' and 05'...HC(8) similar to 
the starting geometry (Fig. 2A) and PCILO preferred conformation 
of m^G obtained without glycosyl torsion angle rotation (5^=16) [18]. 
Hence, this s-tmns conformation of m^G would form canonical 
Watson-Crick base pairing interaction with C44 and non-canonical 
Watson- Crick base pairing with A/U44 in order to provide 
structural stability to the tRNA molecule during protein biosynthesis 
process similarly as observed in earlier conformational and sequence 
analysis studies [18]. 



Stabilization ofs-cis m^G26 conformation 

The MD simulation snapshot structure selected at 2 ns (Fig. 3A), 
and average structures at 5 to II ns (Fig. 3B) and 13 to 19 ns (Fig. 
3C) shows s-cis orientation for methyl substituent of m^G which 
point towards the N(I) atom of guanosine. 

This orientation of N^-methylguanosine allows non Watson- 
Crick base pairing with adenosine (A44) and uracil (U44) instead of 
usual Watson-Crick base pairing with cytosine (C44) at the hinge 
region of tRNA. This s-cis conformational behavior of m^G was also 
noticed in our earlier study [18] and in tRNA^^' crystal structure 
when m^G is present at 1 0th position [34]. The usual Watson-Crick 
base pairing between m^GI0:C25 is not feasible when m^G prefers s- 
cis orientation as observed in crystal conformer (PDB ID: IOB5) 
[34] instead it would form other non Watson-Crick base pairing 
interactions with A44 and U44 at the hinge region of tRNA. 

The geometrical parameters for torsion angles and hydrogen 
bonding interactions analyzed from average and snapshot structures 
are given in table I. The snapshot structure for m^G (Fig. 3 A) taken 
at 2 ns prefers s-cis conformation due to change in a torsion angle 
which deviates from 180° to 29°, while other torsion angles p and % 
diverges to great extent from initial structure and are in close 
agreement with crystal conformer lOBS.pdb [34]. This structure is 
stabilized by the hydrogen bond between N(3)...HC2' (Table I) as 
found in earlier conformational study of m^G [18]. Average structure 
(Fig. 3B) chosen for the period 5 to II ns when a torsion angle 
flipped by 179° as compared to preferred structure of m^G (Fig. 2 A). 
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Table 1. Geometrical parameters for torsion angles and hydrogen bonding interactions for average and snapshot structures after MD 
simulation. 



Modified 
nucleoside 


Average structure 
at time (ns) 


Torsion angle (degree) 


Atoms involved (Atom 1 - 
Atom 2 - Atom 3) 


Distance atom pair 
Atoml- Atom -2(A) 


Angle Atom 1- Atom 
2- Atom-3 (degree) 


Figure 
Ref. 




PCILO most stable 
structure [18] 


a=180°, P=60°, x=286° 


N(3). 
N(3). 


..H-C2' 
..H-C3' 


1.992 
2.269 


117.02 
113.78 


2A 




0-1 


a=180°, |3=333°, x=131° 


N(3).. 


.H-C(IO) 


2.776 


97.41 


2B 




2 


a=29°, |3=254°, x=340° 


IN p j . 


. . ri-v_>Z 


2.503 


125.19 


3A 




3-4.5 


a=180°, P=180°, x=12° 


N(3). 
N(3).. 


..H-C2' 
.H-C(IO) 


2.898 
2.756 


103.42 
93.21 


2C 




05-11 


a=r, P=180°, x=77° 


N(3). 


..H-cr 


2.705 


106.51 


3B 








N(3). 


..H-02' 


2.046 


158.15 






12 


a=l 74°, P=316°, x=333° 


N(3).. 
N(3). 


.H-C(IO) 
..H-C2' 


2.537 
2.902 


91.65 
93.55 


2D 




13-19 


a=355°, P=176°, x=359° 










3C 








05'.. 


.H-C(8) 










20 


a=179°, |3=119°, X=14° 


N(3).. 
N(3). 


.H-C(IO) 
..H-C2' 


2.838 


93.97 


2E 




PCILO most stable 
structure [18] 


a=0°, P=60°, 7=60°, x=286° 


N(3)...H-C2' 
N(3)...H-C3' 


1.992 
2.269 


117.02 
113.78 


5A 




2-3 


a=348°, P=178°, 7=179°, 










5B 




19-20 


a=99°, P=178°, 7=176°, 
X=359° 


05'.. 


.H-C(8) 


1.587 


106.52 


5C 



The obtained average structure maintains distal conformation for 
methyl substituent of guanosine and gets stabilized by N(3)...HCI' 
interaction which was not observed in s-trans conformer of m^G. Last 
average structure (Fig. 3C) was taken within the range of 13 to 19 ns, 
showing s-as conformation of m^G. The s-ci's conformation is 
obtained due to change in a torsion angle from 180° to 355°. This 
average structure (Fig. 3C) shows deviations for torsion angle p by 
116° and X t>y 73°. Obtained average structure (Fig. 3C) shows 
similar values for torsion angle a and % as compared with crystal 
structure IOB5.pdb [34]. A large deviation around the torsion angle 
P is due to fluctuations from s-trans to s-cis conformation by rotating 
C-N bond of methyl group. 

Fluctuations in torsion angles of JSF -methylguanosine (m^G) 
during MD simulation 

Analyses were also made for torsion angles and hydrogen bonding 
interactions of m^G during 20 ns simulation period (Fig. 4). The 
torsion angle a fluctuates periodically between s-trans (+180°) and s- 
cis (0°) rotamers of m^G during total simulation period (Fig. 4A and 
Table I). For simulation time 0 to I ns, 2.8 to 4.7 ns, 1 1.5 to 12 ns 
and 20 ns (Fig. 4A) torsion angle a prefers s-trans orientation which 
is supported by weak interaction between N(3). . .HC2' (Fig. 4D) 
and N(3). . .HC(IO) (Table I). Orientation of a torsion angle favors 
the usual Watson-Crick base pairing of m^G26 with C44 and unusual 



non- Watson-Crick base pairing with A/U44. Whereas, during 
simulation period I to 2.8 ns, 4.8 to IL2 ns and 12 to 19.8 ns 
methyl substituent of guanosine prefers s-cis orientation which is 
stabilized by N(3)...HC2' and N(3)...HCr hydrogen bonding 
interactions. 

The s-cis conformation of N^-methyl substituent has also been 
observed in crystal structure when m^G present at 1 0th position in 
tRNA^^' [^4]* This orientation allows the non Watson-Crick base 
pairing between s-cis m^G with A/U44. Similar results were found in 
earlier conformational energy calculations performed over m^G [18]. 

Torsion angle p (Fig. 4B) maintains starting geometry +180° 
[18] with small fluctuations at dl 60° as found in crystal structure 
6TNA.pdb [9], lEHZpdb [33] and IOB5.pdb [34]. The glycosyl 
torsion angle (x) (Fig. 4C) fluctuates at ±30°, ±120°, ±180° and 
favors the respective anti{\!^ to 9.2 ns, IL2 to 18.9 ns and 19.2 to 
20 ns) and syn (0 to L7 ns, 9.3 to 10.2, 10.8 to II ns and I9ns) 
conformation for m^G. The m^G is preferably stable at both syn and 
anti conformation, which allows usual (Watson- Crick) as well as 
unusual (non Watson-Crick) base pairing with C/A/U44. The syn 
conformation of s-trans and s-cis m^G is supported by 05'...HC(8) 
hydrogen bonding interaction whereas in and conformation of s-trans 
and s-cis m^G, is held by N(3)...HC2' and N(3). . .HC(IO) (Table 
I) during MD simulation. Hydrogen bonding between 05'...HC(8) 
is varied in accordance with the fluctuations found in glycosyl torsion 
angle {y^ of m^G26 during simulation period. 
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Figure 4. Molecular dynamics (MD) result: (A) Showing fluctuations in a torsion angle. (B) Fluctuations in (3 torsion angle. (C) Fluctuations in x torsion angle. (D) 
Fluctuations in hydrogen bonding between N(3)-HC2' (E) 05'-HC(8) 



Molecular dynamics (MD) simulation study oflSF-lSF dimethyl 
guanosine (m^iG) 

In order to see solvation effect on N^-N^ dimethylguanosine 
explicit molecular dynamics simulation study of 20 ns has been 
performed over the PCILO preferred conformation (Fig. 5 A) [18]. 
To confirm the conformational behavior of m^zG we have analyzed 
two different average structures taken at 2-3 ns (Fig. 5B) and last 
1000 ps from 19-20 ns (Fig. 5C), their geometrical parameters are 
listed in table I. 

The average structure obtained from 2 to 3 ns (Fig. 5B) prefers 
distal conformation for m^zG and prevent Watson-Crick base pairing 
with C, instead it would prefer non-canonical Watson-Crick 
interactions to pair with A/U44. Compared with crystal conformer 
(lEHZ.pdb) [33], average structure retains quite similar torsion angle 



values for a, p and The average structure selected at last 1000 ps 
(19 to 20 ns) does not show much difference as compared with 
earlier average structure. The only difference between these two 
average structures is variation around a torsion angle which is 
positioned to 99° (Fig. 5C) from its preferred value (Fig. 5 A). Due to 
this small change in conformational property of m^zG, dimethylation 
of guanosine has almost no influence on pairing with either A or U, 
because the position of guanosine has no impact on these base 
pairing interactions. The m^zG-A pair can be formed with little 
hindrance, because, even though the methyl groups of m^2G26 are in 
the plane of the base, as in case of yeast tRNA^^^ they take part in a 
propeller-type arrangement with the base. The m^zG-U pair would 
not be affected by any conformational arrangement of the methyl 
groups similarly as discussed in [18]. 
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Figure 6. Molecular dynamics (MD) result: (A) Showing fluctuations in a torsion angle. (B) Fluctuations in (3 torsion angle. (C) Fluctuations in y torsion angle. (D) 
Fluctuations in x torsion angle. (E) Fluctuation in hydrogen bonding between N(3)-HC2' (F) 05'-HC(8) 



Fluctuations in torsion angles of ISPy ISP- dimethylguanosine 
(m^2G) during MD simulation 

The a torsion angle fluctuates periodically in between +180° or 
d: 60° (Fig. 6 A and Table I) over the 20 ns molecular dynamics 
simulation period suggesting free rotation around C-N bond. 

Torsion angles p (Fig. 6B) and y (Fig. 6C) retain preferred values 
lb 1 80°, with small fluctuations at +60° throughout the simulation 
period. Glycosyl torsion angle {j^ adopts and conformation during 
the simulation study. Such type of anti conformation for N^,N^- 



dimethylguanosine was confirmed through crystal structure 
(lEHZ.pdb, lEVV.pdb, lOBS.pdb 6TNA.pdb). 

Conclusion 

The regular periodical fluctuations around the bond C(2)-N(2) 
of m^G was observed throughout the 20 ns molecular dynamics 
simulation, which confirms the existence of iso-energetic s-cis or s- 
tmns rotamers of m^G. These iso-energetic rotamers interconvert 
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easily during the simulation period. These results are in favor with 
preferred and alternative conformations of m^G obtained by our 
earlier conformational energy calculations [18] as well as crystal 
structure (lEHZ.pdb [33] and lOBS.pdb [34]). The periodical 
fluctuations of s-tmns to s-cis and vice versa could be possible due to 
free rotations around the C-N bond of methyl group. According to 
tRNA sequence analysis [18] and this MD simulation results we 
would like to say that m^G26 can form three different canonical as 
well as non-canonical Watson-Crick base pairing interactions with 
other bases. Such base pairing may be summarized as i) an usual 
Watson-Crick base pairing of m^G26-C44 where the methyl 
substituent must be in s-trans orientation, ii) non Watson-Crick base 
pairing between m^G26-A where the methyl substituent is likely to be 
s-cis orientation, and iii) non Watson-Crick m^G26-U base pairing 
where the methyl group can adopt one of them, i.e. s-cis or s-tmns 
conformation. These results reveal that m^G is equally stable as either 
the s-cis or s-tmns rotamers and the rotational preference of methyl 
group may be specific to the sequence context reliant upon which face 
of the base contributes in hydrogen bonding. Thus, MD simulation 
results confirm that the N^-methyl group of m^G26 may prefer 
energetically two stable rotamers, i.e., s-tmns m^G26:C/A/U44 and 
5-a5 m^G2 6: A/U44 as found in earlier results [18]. 

Similarly, the presence of two methyl groups unlike in case of 
single methyl in m^G virtually eliminates the possibility of pairing 
with C and, indeed, m^2G26 pairs exclusively with A or U at position 
44 and is flanked by C27:G43 on one side and the m^GI0-C25-G45 
triple on the other [35]. Hence, these results suggest that the modified 
nucleosides m^G26 and m^2G26 play an important role in tRNA 
folding and may prevent tRNA from adopting the unusual 
mitochondrial like conformation. 
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